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Abstract 

Magnetization ground states are studied in toroidal nanomagnets. The energetics associated 
to the ferromagnetic, vortex and onion-like configurations are explicitly computed. The analysis 
reveals that the vortex appears to be the most prominent of such states, minimizing total energy 
in every torus with internal radius r > 10 nm (for Permalloy). For r < 10 nm the vortex remains 
the most favorable pattern whenever R/lex ^ 1-5 {R is the torus external radius and lex is the 
exchange length), being substituted by the ferromagnetic state whenever R/le.x ^ 1-5. 

1 Introduction and Motivation 

Nanomagnetism has become an important area for both basic and applied research in the last 
years. Such an importance is greatly due to recent advances in the synthesis and characterization of 
magnetic samples of very small sizes, comprising micro and sub-micrometer scales. Nanomagnets of 
several shapes and sizes have appeared and have been intensively investigated, namely for potential 
applications in logic memory, data storage and highly sensitive sensors [T]. Such an interest has been 
renewed and reinforced by virtue of a very recent proposal of using biofunctionalized magnetic- vortex 
microdiscs for triggering cancer therapy [2]. Often, such applications are based upon new magneto- 
electronic mechanisms which preclude on demand control of the magnetization stable configurations, 
like single-domains, domain-walls, vortices, etc. Among other, cylindrical and rectangular-type shape 
samples appear to be the most studied ones, once their relatively simple geometries make easier the 
fabrication and eventually the comparison between experimental and theoretical findings. For in- 
stance, in a thin enough sub- micron sized disc of soft material, like Co or Permalloy (FeNi), the 
remanent (at zero applied field) ground-state is observed to be a single- vortex of rotating dipoles with 
a central core where they revolve against the disc plane for regularizing exchange cost [3^ . As the disc 
thickness is raised the vortex is progressively distorted being eventually substituted by a single-domain 
configuration [?]. In turn, thin rectangular- like samples generally present Landau patterns as their 
remanent ground-states [S]. 

It is widely recognized that the total energy of a magnet with an arbitrary magnetization pattern 
is not easy to be computed, even assuming simple geometries and negligible anisotropics. With these 
simplifications, exchange and magnetostatic contributions remain to be evaluated [HIT]. The latter 
one is generally much more complicated once it comes from long-range dipole-dipole interactions. 
For special configurations, like vortices and some aligned states, and/or simple geometries, like those 
curvatureless, say, planar-type and cylinder, its evaluation is greatly simplified and may be sometimes 
exactly computed. For more general situations, the available tools for doing that often relies on ap- 
proximate methods, e.g., assuming the nanomagnet as a continuum sample (analytical treatments) 
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and/or performing numerical evaluation of the energy [5]. For the specific cases of the square- type and 
cylindrical shapes, considerable amount of information are available with relatively good agreement 
between experimental and theoretical studies [8]. For other geometries, a detailed knowledge is still 
lacking, although some works have dealt with nanosized magnetic spheres [9], wires [1^, tubes [11] 
and rings [12l|13l[I3. 

Whenever designing novel magneto-electronic mechanisms based upon vortex-like states, namely 
those for logic memory and/or data storage, an issue of prime importance concerns miniaturization 
of each nanomagnet and consequently of the whole array comprising a large number of them. Since 
current transistor-based memory elements have reached below 100 nm per element, the departure to 
nanomagnetic vortex-based mechanisms should offer real potentiality of supporting such stable states 
at much smaller scales, say, at least around 20 — 30 nm, or less. There lies a trouble with soft material 
nanodiscs, where vortex stability is lost if the disc radius gets lower than ^ 45 nm (value for Permalloy 
disks with thickness 10 nm) |l5j . To partially remedy the situation, one may introduce a (centralized) 
hole into the disc, obtaining a nanoring, where a vortex-state may be stabilized at smaller radius 
20 nm, for internal radius and thickness about 10 nm) [13]. Physically, the presence of the hole 
becomes the core formation no longer necessary, diminishing exchange cost considerably and giving 
extra stability to vortex configuration even at smaller scales, see Refs. [TCKT^. However, disappear- 
ing with the core, we immediately loose the vortex polarization, remaining only its chirality (sense 
of dipole rotating, clock or counter-clockwise) as bit storage media. Additionally, dipole dynamics 
at the edge borders may still jeopardize vortex stability due to large fluctuations of surface magnetic 
charges, (Tmag = M ■ h, once at these places the normal vector n is not well-deflned (in other words, 
curvature blows up at disc and hole edge borders) . Such fluctuations may be somewhat bypassed and 
their effects minimized if the centralized hole is large so that the nanoring width is sufficiently small, 
less than a domain- wall size, say, around a couple of exchange length (£ex ~ 5.3 nm, for Permalloy). 
Doing that we also ensure that no vortex core can be formed in the ring [14] . 

Here, we would like to consider a number of magnetization patterns lying in the toroidal manifold. 
Specifically, we compute the net energy (exchange -I- magnetostatic) of ferromagnetic (single-domain, 
SD), vortex and onion- like states. Our main results may be divided according nanotori internal radius, 
r: for small values, r < 10 nm, SD appears to be the ground-state if R/iex ^ 1-5, being substituted 
by a vortex configuration whenever R/£ex ^ 1-5. Thus, even in a very small torus, say, with r = 1 nm, 
we may have a stable vortex if i? ^ 7.5 nm. On the other hand, if r > 10 nm only the vortex pat- 
tern appears to minimize the net energy, despite the aspect ratio. Such a vortex stability, even at 
very small scales, should be attributed to the toroidal geometry, where curvature varies smoothly, 
preventing large fluctuations of surface charges at the edges, as may occur in nanorings. A direct 
consequence concerns the possibility of great improvement in miniaturization using stable vortices 
lying in tori as media for binary storage and logic. Note that such vortices are coreless, leaving only 
their chirality at our disposal. However, the vortex core absence is an important property when we 
consider applications where non-interacting magnetic elements concern. For example, the reduced 
dipole-dipole interaction may, in this situation, increases the vortex-state stability for small intertori 
separation distances, whenever compared to its interdiscs counterpart. Such a feature may be impor- 
tant for fabricating arrays with higher density, using nanotori as the basic elements. For applications 
in cancer therapy, the vanishing magnetization in remanence is significant because this reduces the 
long-range magnetostatic forces responsible for particle agglomeration. 

To our best knowledge, there is so far no work reporting the fabrication of magnetic nanotori. 
The best we have found in the literature concerns some possibilities of producing such structures 
using their carbon-made analogues as a sort of substrate where a suitable magnetic material could 
be formed either by capping the surface with a thin layer or by injecting the material inside to form 
a nanotorus encapsulated by the layer 181 . By means of this latter technique, it has been reported 
in Ref. |19j Permalloy-made nanowires with (internal) radius around 10 nm (other possibilities for 
nanowires fabrication may be found in Refs. |20j). Once experimental techniques in nanomagnetism 
have developed rapidly it is expected that, among other, magnetic nanotori may appear in the near 
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future. 

At some extent, this work represents a further step from that presented in Ref. |21| . where we 
have studied sohtons and vortices in the toroidal geometry considering only the exchange energy. Here, 
for dealing with nanostructures their magnetostatic contribution must be also taken into account. We 
outline our article as follows: in Section 2, we describe the theoretical model, the torus geometry and 
the magnetization states under consideration. Section 3 is devoted to the analytical results concern- 
ing the energy of these magnetic states. In Section 4, we evaluate numerically the previous findings, 
followed by a suitable discussion about them. Finally, we conclude our work and present a number of 
prospects for future investigation. 



2 The model and the toroidal magnetization states 

In our approach, a magnetic nanotorus is assumed to be a continuum material where the magne- 
tization vector has the same magnitude throughout it, but can vary in direction from point to point. 
Its energy comes from short-range exchange (E'ox) plus long-range dipole-dipole (magnetostatic, E'mag) 
interactions, once we assume vanishing anisotropics, i.e., we are explicitly dealing with isotropic ma- 
terials, like Permalloy (FeNi) and Co. Therefore, the ground-state will be the configuration which 
minimizes -Enet = Ee^ + i?mag , which generally depends upon the torus geometry. Explicitly, we have 
that the exchange term is given by the continuum version of the Heisenberg model over the whole 
volume [2T| : 




where A accounts for the interaction strength among nearest- neighbor spins (dipoles). Once we are 
considering a ferromagnetic sample, ^ > 0, so that -Bex > 0; rn{r) — {mx,my,mz) — M/Ms is the 
normalized magnetization and Ms is its satured value. For Permalloy, one has A « 1.3 x 10~^^ J/m 
and Ms « 8.6 x lO^A/m, so that, £cx — \/^A/JloMJ « 5.3 nm. Such a length accounts for the 
relative importance of each contribution, in such a way that exchange dominates at small distances, 
d < ^ex, whereas magnetostatic takes over at larger ones, d > £cx- This latter energy comes from the 
interaction experienced by a given dipole due to the remaining ones throughout the magnet. With 
the additional assumption that no electric current is present, J = 0, we formally obtain: 

Em^g = ^ m(f) • V^f)dV ^ 0, (2) 

with the magnetotatic potential given by: 

where Pm(f^) — V' • rn{r') and (Tm{r') — n' ■ m{r') are the volumetric and superficial densities of 
effective magnetic charges, respectively. 

Since our interest concerns ferromagnetic tori, it is instructive to give a survey of its shape with 
the main features. Topologically, an ordinary torus is the simplest two-dimensional closed manifold 
having genus- 1, a single hole, like depicted in Fig. [TJ Depending on the relative sizes of its internal, 
r, and external, R, radii, three distinct standard types are obtained: ring, horn, or self-intersecting 
spindle torus, according R>r,R — r, or R<r, respectively (see, for instance, Ref.[lT] and related 
references therein). The commonest one is the ring torus, resembling a donut, which will be the 
type studied here (the other ones, although mathematically interesting, seem to be rare in Nature). 
Geometrically, a torus is a surface whose curvature smoothly varies along the polar angle, 0, from 
negative to positive values. For describing its geometry some coordinate systems are frequently used, 
among them, the spherical-type relates the Cartesian to the toroidal angles, as below: 
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Figure 1: [Color online] Commonest kind of torus, the ring torus, and an useful eoordinate system 
over it, {0,(j)) G [0,27r]. Parameters r and R are the internal and external torus radii, respectively, 
whereas R/r is the aspect ratio. 



with {0,(p) E [0, 27r] being the polar and azimuthal angles, respectively (see Fig. [T]). Another useful 
coordinate system, along with a number of additional relations used for evaluating some quantities 
obtained in Section 3, may be found in Appendix A. 

Let us now discuss the main aspects concerning the relevant magnetization states whose ener- 
getics will be evaluated in Section 3. First, let us recall that the ferromagnetic state is reached when 
all the dipoles align along the same direction, leading to magnetization saturation. By virtue of the 
parallelism between neighbor dipoles, their exchange energy vanishes identically. In addition, no vol- 
umetric magnetic charges take place for this pattern, pm = V • m = 0, so that, only surface charges, 
Cm = n ■ m ^ 0, contributes to its magnetostatic energy. Therefore, the results concerning a toroidal 
surface may be straightforward applied to a volumetric torus. In section 3, we shall compute the total 
energy for such a pattern lying along z and xy-jAane, given respectively by, mjr^ — z and rnj^^y = x 
(any other direction along xy is energetically equivalent). As a matter of fact, we shall see that in any 
torus the mjr^^-state is energetically favorable than its z counterpart. 

Another relevant state to be considered is that which minimizes magnetostatic term, which is 
achieved by vanishing both pm and Um- This corresponds to a curly pattern around the torus genus 
such that its spatial distribution ensures V • m = 0, besides being perpendicular to the surface at 
every point, so that n • m = 0. Such a vortex-state may be described by m„ = ±0, where is the 
unit vector along the azimuthal angle, whereas ± is associated to the vortex chirality, counter- or 
clockwise, respectively. Furthermore, it is easy to obtain that \J{nix + my) leading to a non-zero 
exchange contribution, given by eq. (P), to be evaluated in Section 3. 

Additionally, we also consider a more complex structure, known as an onion-like state. It may 
be described by a magnetization depending on the azimuthal angle, rfiQ — mr{(f>)R + m^{(f>)(l), where 
R = X cos (j) -\- y sin and \22\ : 




X 



(i? + r sin0) cos 0, y = (i? + rsin^) sini/) and 



z = r cos 9 , 



(4) 





< < 7r/2 
tt/2 < (f) < n 
TT <(/)< 37r/2 
37r/2 < (/«< 27r 



with f{q, (p) = /((/)) = cos''{(j>) E [—1, 1]. The real parameter q> 1 accounts for the relative deviation 



Figure 2: [Color online] Some configurations considered in this work. From left to right, one can see 
the single-domain (along z), the vortex, ferromagnetic and onion (with q — 7; see text for details) 
patterns, respectively (the last three ones with the dipoles lying along xy plane). 



of the pattern compared to the ferromagnetic state, "nijr^^ = i, formally recovered as long as 9 — )■ 1 
(Fig. [2] shows a number of such configurations). Thus, as higher is the q value more is the respective 
onion pattern distorted compared to the ferromagnetic state. These are the relevant configurations 
usually studied as possible ground-states in other nanomagnets, like discs, rings and spheres. Our 
intention here is to proceed with a similar plan in order to determine which of such states minimize 
the net energy for a given ring torus, whose specific geometry is dictated by its aspect ratio, R/r. 



3 The energetics of the configurations: analytical results 

3.1 Exchange term 

As already said, ferromagnetic pattern has vanishing exchange cost, once magnetization does not 
vary throughout the sample, so that Vm = 0. 



Vortex exchange energy may be analytically computed within our approach based upon the 
continuum version of the Heisenberg model, eq. (HJ, giving: 



y/b^ +t^- \b\ 



(6) 



where P — E? 



> and \t\ = Vr 



e (0, r] and c is the hole radius of the hollow torus 



{c — represents the bulk torus), 
counterpart, £^ 



Expression above may be also obtained from its two-dimensional 
A-K^Jr/y/R^ — r^, by integrating over r from c to r, see Ref. [21] [A = kJ/a, where 
a is the lattice spacing parameter, J is the 2D coupling between the neighbour spins, the stiffness 
constant with energy dimension, and k a real parameter depending upon the lattice structure, e.g., 
k — 1 for simple cubic and k = 2 for body-centered cubic crystals). 



Analogously, the onion-state exchange cost may be evaluated to give: 

Eg^E^,{liq)-l), (7) 

where 

2 r^'^ 1 

with I{q) > 1. Indeed, X{q) = 1 for g = 1, so that the onion-like pattern recovers the ferromagnetic 
state along the xy-plane. In addition, I{q) > 2 whenever q > 7, so that such configurations (with 
q > 7) cannot occur as ground-states, once their total energy is higher than the vortex one, as one 
can realize from Fig. O 



df{q,<P) 
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Figure 3: [Color online] Exchange energy as the function of the torus aspect ratio, R/r, for a number 
of possible configurations. Namely, note that whenever q > 7 the onion-type is more energetic than 
the vortex state. 



3.2 Magnetostatic term 

Once no effective magnetic charges appear in the vortex state, its total energy comes solely from 
exchange cost, computed above. The remaining patterns have magnetic charges lying in the torus 
surface and/or volume, making their energetics rather involved. Below, we shall present only the re- 
sults and their main features, leaving the details for the respective Appendices. In every case, we start 
by determining the magnetostatic potential, from what the magnetic field emerges hy H — — V<i>„i, 
as usual. Our procedure is based upon the evaluation of the 3D Green-function in terms of toroidal 
harmonics, involving Legendre functions of half-integer orders, P'^±i/2 ^^^^ Qn±i/2^ of first and second 
kinds, respectively (we refer reader to Appendix B for details). 

Let us start by the simplest of such cases, the ferromagnetic (or single-domain, SD) along z-axis, 
rfijr^ = z, in which only a^n appears. We formally obtain (see Appendix C.l for details): 

EL, = ^^^l^f^ sinh^ ao £ n^P^^i/^g^-i/^ {qU/,) ' , (8) 

n— 1 

where we have omitted the superscripts fc = in the Legendre functions and also defined P^(cosh ao) = 
and Q'^{coshao) = Q'^, where coshao — R/r is the torus aspect ratio. Although we have explicitly 
assumed a bulk torus, the result above readily gives its hollow counterpart (we shall return to this 
point later). 

Below, we present the magnetostatic contribution associated to onion-like patterns, whose special 
case 9=1 represents the single-domain along the plane, SDa;^. Now, both surface and volumetric 
charges take place, making the determination of very intricate. Even for this involved case, we 
have obtained a formal analytical expression for every q, as below: 



fe=0 n=0 \ 2' 



(9) 
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with the definitions: 

Jl^ Jfe((^')"^r(^')V 

Jo 

<3«-i/2(cosha) 



T" — 
>-'fc — 



2tt 



sinha 



da, 



Above, we have used the coordinates (a, (3, if), presented in Appendix A (namely, we have that 
R = 6tanhQ;, r = fosinha, and (j) = f), whereas Am Jk{f'), B^ and, Gniv') are given in Appendix 
C.2. Further evaluation of energy (jO)) demands numerical integration with (i?, r) as inputs. 



By taking g = 1 in Eq. ([9]) we obtain the magnetostatic energy associated to SDa;y, as below 



TpJ~'xy — 

mag 



E 

n=0 



where we have used the property —i — ItQ,\-i/2 



tri-^n^n-l/2^n-l/2' 

Qn-i/2 (see Ref. [23", p. 959). 



(10) 



Let us return to the hollow torus, with thickness t. The potential $m may be readily computed 
through linear superposition of those associated to two concentric tori with radii r — t and r, having 
opposite magnetization on their surfaces. After that, an analogous procedure yields the magnetostatic 
energy. To make use of previous results, let us take coordinates (a, /?, (p), so that two arbitrary tori, A 
and B, are parametrized by (a^, /5, v) and {aB,f^', (f'), respectively. Then, we have cosh(a^) — R/ta 
and cosh(aB) = R/rg, and by assuming r^i > r^, we obtain, for the ferromagnetic state along z: 



hollow 



128^0^1 
9 



OO 

11=1 



h\ Sinh^ aAPtl/2Qi-l/2 {QI-I/2 



-b% sinh^ aBPti/2Qn-i/2 {Q'n-1/2) 
whereas for that along the plane, we have: 



(11) 



mag 



32moM|^ 



hollo 



n=0 



l3 plA r)-^^ 

"A^n-l/2^n-l/2 



cosh a A 



/2 



sinh a A 



l3 pis r^~lB 



cosh a B 



"^n+1/2 

sinh a b 



\1B 



^\-l/2 

H 2 N 



Q 



2A 

n-1/2 

sinhck^ 



-1/2 



/2 



i?2 



sinh a b 



, while Pj^ {cosh a AiB)) = Pj^^'-^^ and Qi^{cos\iaAiB 



(12) 



Above, we have defined bA.B - 

Qi^^(^\ The second terms in these expressions, proportional to &b, identically vanish as long as 
tb — ^ and ta = r, as expected. The same may be done for the onion-like state, but due to its 
length, we shall not quote it here. 



4 Numerical results and discussion 

Although we have obtained analytical expressions for the magnetostatic energy of the desired 
magnetization in the torus, Eqs. (|5))- (ll0p . an immediate analysis of their physical content is not 
clear due to their complicated forms. Therefore, to extract relevant results we should evaluate them 
numerically. For doing that, we employ the Fortran subroutine code DTORHl [24j, which gives the 
values of Pl^_^^^(cosha) and Q'^^^^^i^osha) for a given torus with radii R and r. We also adopt the 
values of A and Ms for Permalloy, already given in Section 2. 
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Figure 4: [Color online] The magnetostatic energy of the ferromagnetic states along z and xy against 
tori volume {y = 2TT^r^R). Above, we have fixed r — 5nm whereas R varies from 5.05 nm to 2500 
nm. Clearly, SDxy is energetically favorable compared to its counterpart along z. 



First, we compare the energies associated to the ferromagnetic (single-domain) states along z and 
xy, which may be summarized by the expressions below (V = 2Tr^r^R is the torus volume): 



^mag ^ mag ^ ^ on 

~ 8 ' MoA^I ^ 4 ' ^ ' 

clearly showing that SD^y costs nearly half of the energy from its counterpart along z (Fig. |4|). Thus, 
hereon only this ferromagnetic pattern need to be considered. Next, recalling that it is a particular 
case of the onion-like state, with q = 1, it is important to determine the values of this parameter that 
minimizes total energy for a given torus. This may be formally achieved by solving dE/dq = 0, and 
the results are summarized in Fig. [SI where we present the net energy of vortex, onion and SDxy 
states against R/r. For each graphic, the onion state loses its stability at the point O, once it is more 
energetic than the vortex state for larger R/r. At this point, numerical results show that q « 1.7 
minimizes the total energy in a torus with i? « 8 nm and r — 1 nm (for larger R, the ground-state 
appears to be the vortex pattern) . We can also note that the total energy of the onion state, at the 
point O, is only about 3.2% less than its ferromagnetic (SD^jj,) counterpart (indicated by the point 
F in Fig. [5]). This fact must be associated with the little deviation from parallel direction of the 
magnetic dipoles for q = 1.7, as we can see in Fig. [6] As r is raised, such an energy difference between 
those configurations is diminished, as illustrated in Fig. [5] (on the right) where r — 2 nm has been 
taken. Therefore, for saving time and computational efforts, we shall adopt g = 1 as a very good 
approximation for evaluating the net energy, say, for our purposes it is suffice to compare vortex and 
SDxy states energetics. 

To complete our analysis, we should compare vortex and SDxy energetics to decide what is the 
magnetization ground-state for a given torus. Figure [7] summarizes our main findings, where we have 
indicated the energetically favorable state for a given torus aspect ratio, R/r. Considering only the 
physically suitable region, with R> r, we may realize that the vortex appears to be the ground-state 
in every torus with r > 10 nm [p = r/t^x ^2). In those tori with r < 10 nm vortex still minimizes the 
total energy whenever R > 10 nm (^ = R/£ex ^ 2) whereas SDxy is energetically favorable for smaller 
R. Mainly in the region p < 1 and ^ < 1.5 the true ground-state is not SDxy {q = 1) exactly, but 
slightly different onion-like patterns with 1 < q < 1.7, whose energies are, at best, 3.2% below that 
the q — 1 state. Therefore, the transition curve presents a confidence around this percentage in that 
region. 
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Figure 5: [Color online] The total energy corresponding to vortex, onion and SD^y (9 = 1) for r = 1 nm 
(left) and r = 2nni (right). In each case, the point O marks the critical R/r value above which the 
onion state is no longer stable, say, its energy is higher than the vortex-state {F is the point at which 
vortex and SDxy have the same net energy). The highest difference between the energies of the onion- 
like (with q ^ 1.7) and SD^jj, states occurs for a torus with r — 1 nm and R ~ 8nm, and reads « 3% 
(see discussion in text). 




Figure 6: [Color online] Illustration of onion-like configurations with q — I (ferromagnetic pattern; 
left) and q = 1.7 (right). As also indicated by their respective energetics, these states are very similar 
each other (see text for further details). 
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Figure 7: [Color online] State diagram for nanotori = R/tex and p — r/if-x)- For every tori with 
r > 10 nm (p > 2) the vortex pattern appears to be the ground-state. At smaler radii, ^ 10 nm, 
vortex remains the ground-state in tori with R > 8nm > 1-5), being substituted by SHxy at smaller 
R. The transition curve is precise within sa 3% of confidence (for p "^1 and ^ < 1.5). 



It is now instructive to compare our state diagram with its ring counterpart, for instance, that 
presented in Ref. [13 (unfortunately, in this work authors do not provide the energy values for a better 
comparison among their and our results). Among other, a remarkable contrast concerns how small 
these geometries may support a vortex-like as the ground-state, which may occurs even in very small 
torS say, with r sa 1 nm and i? > 8 nm, whereas in a nanoring at least an external radius, > 10 nm, 
with a central hole and thickness no less than 8 nm and 2 nm is demanded. More specifically, if 
r = 1 nm (thickness 2 nm), vortex pattern is stable for nanotori with volume V ~ 150 nm'^, whereas 
for nanorings with the same thickness, this stability demands a volume not less than V w 225 nm'^ 
(around 50% higher than its nanotorus counterpart). Therefore, by virtue of its smooth curvature, 
the toroidal geometry offers a remarkable physical support to sustain miniaturization along with a 
single-vortex magnetization ground-state. In addition, thinking about an array of such nanotori, 
a system of practically non-interacting elements can be achieved provided that a separation > tf,x 
among neighbor nanotori is ensured. Although such a non-interacting feature may be obtained with 
nanorings, the nanotori has the advantage of minimizing considerably the magnetostatic fluctuations 
due to magnetic surface charges at the edges. 

Now, we discuss upon the hollow torus case, whose vortex and single-domain energies are given 
by Eqs. ^ and (|12|) . respectively; two relevant state-diagrams are shown in Fig. [H Namely, note 
that a higher value of R is demanded for ensuring the vortex ground-state than in a bulk torus with 
the same r. Thus, in a hollow torus SD^^y pattern appears to be more favorable as the torus thickness 
is diminished. However, at least for the cases considered, say, a torus with a fixed thickness t = 1 nm 
(FigEl on the left), and with varying thickness according t — r/2 (FiglH on the right), they are 
qualitatively similar to the bulk torus, with small quantitative differences. 



^We assume that our classical analysis remains approximately valid at such scales, of a few atoms along the internal 
diameter, where certainly quantum effects could be important. If we consider r > ^cx ~ 5 — 6nm, our statement about 
vortex stabilization is more reasonable within our approach. 
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Figure 8: [Color online] State diagrams for the hollow torus. On the Left, a fixed thickness t = 1 nm 
and r G [2, 20] nm are taken; on the Right, the thickness varies as t = r/2. Namely, note that vortex 
stability tends to be (slightly) diminished as thickness is lowered, say, a thin hollow torus favors 
single-domain as its ground-state. 

5 Conclusions and Prospects 

We have studied the energetics of a number of magnetization configurations lying in the geometry 
of a toroidal ferromagnet at nanoscales. From our analysis vortex-like emerges as the most promi- 
nent of these patterns, being observed to remain stable even in very small tori, with r « 1 nm and 
R > 8nm. Actually, onion-like configurations, with 1 < (7 < 1.7, which are practically equivalent to 
ST)xy, takes place only for p = r/iex ^ 1 and ^ = R/iex ^ 1-6. This is an important result as long 
as miniaturization of basic magneto-electronic elements comprising a single- vortex concerns, once this 
configuration is stable in nanotori much smaller than in their ring and disc counterparts. 

Under investigation are some possible mechanisms to induce switching between the two chirality 
states of such a vortex and, perhaps more interesting, the study of an array of such nanotori com- 
prising a vortex in each one. In this situation, if the tori are spaced more than an exchange length 
from each other, we have a system of practically non-interacting elements. Indeed, even possible in- 
teractions induced by external agents (applied field, etc), which may deform the vortex, are expected 
to be minimized due to the toroidal smooth geometry. Additionally, such a smoothness is expected 
to prevent large magnetostatic fluctuations coming from surface magnetic charges, as may occurs in 
discs and rings straight corners. Furthermore, that weak self-interaction characteristic among the ba- 
sic elements could make magnetic-vortex nanotori assembles as interesting systems in future targeted 
cancer-cell therapy 
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Appendices 



A Toroidal coordinate systems 

Besides the spherical-like coordinate system, presented in Section 2, Eqs. the so-called toroidal 
coordinates are defined by means of the following parametric relations: 

sinh a cos sinh a sin sin/3 

u p' y = ° ^' ^' 

cosh a — cos p cosh a — cos p cosh a — cos p 

where 6 is a constant giving the radius of a circle at z = described by a — >■ cxi, say, when 
a — > oo, we have x = bcos(j), y — bsuKp and z = 0; their ranges are ^ a < oo, ^ /3 ^ 27r and 
^ (p ^ 2tt. On the torus surface we have: 

b = y/(R + a)(R — a), coshao = — , (15) 

r 

so that b and ao may be viewed as the geometric radius and the eccentric angle, respectively. 



B Laplace Equation in Toroidal Coordinates; toroidal 
Green Functions 

The coordinate system (|14l) is very useful because it has a known set of solution for the Laplace 
equation, where the method of separation of variables is valid (see Ref. '55' p. 1302): 



$(a,/3, V?) = ^cosha - cos l3A{a)B{l3)C{ip). 
where the functions A, B and C satisfy {C, — cosh a, and n, k constants): 



dA 



da^ 



da 



n 

4 



A = 0, 



(16) 
(17) 
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Figure 9: [Color online] The behavior of some toroidal functions against R/r. 



^+n^B = (18) 

'^ + ec = 0. (19) 
The functions i?(/3) and C{ip) are linear combinations of trigonometric functions, like below: 

Bn{P) = fn cosn/3 + gn s'mnfi, (20) 

Cki^) = Uk cos kip + Wksmkip, (21) 

where /„, gn, Uk and Wk are constant coefficients to be determined from the boundary conditions. 
Since /3 and ip are variables ranging from to 27r, the solutions must be periodic. This condition 
implies that k — 0,1, 2, 3, ... and n — 0, 1, 2, 3, that is, invariance under rotations imposes k 
and n to be integers. 

The function A{a) are the Legendre functions of half-integer order of first kind ^'^^ second 

kind <5^_]^/2' ^^^"^ called toroidal functions or toroidal harmonics, whose asymptotic behavior are 
given by (see Fig.[9|): 



lima^oo 1/2 (cosh a) = cg 
lima^o Qj;_i/2(coshQ;) = 00. 



(22) 



The expansion l/|r — r'\ in terms of Green Function, according to the coordinate system (|14p . is 
given by (See Ref. [25], p. 1304): 



J (cosh a - cos /3) (cosh a' - cosh P') yr-^ yr-^ , 

— —^^{-iftnek cos n{l3 - /3') cos k{ip~ if') 



fc=0 n=Q 
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r [n + k I 2 



ff^n-i/2(cosha<)g^_i/2(cosha>), (23) 



where = (2 — Snfl), Ck = {2 — (5m, o) and r is parameterized in terms of (a', /?', The fact of 
P^_j^^^{cosh.a) and (5^„i/2('^os^'^) ^"^^ related to < and >, respectively, reflects the asymptotic 
behavior of these functions (see Fig. [3]). 



C Evaluation of the magnet ost at ic energy: expressions ([8]), 
dSD and (UnD 



C.l Evaluation of Eq. dS]) 

In the case of the SD^, we have that M — Msz, what gives: 

sinh a sin B 

M-h = Ms— ^, 24 

cosh a — cos p 

so that only k — component contributes to the sum, and the integration in ip is easily evaluated 
to give 2tt. Thus, the substitution of p3|) in Eq. ([3]) yields (fc = superscripts has been omitted): 

6Mssinh^a f^"" sin (3d(3 , — 

= / — TTjV cosh a' -cos /?'>£„ cos n/3-/3) 

27r Jq (cosha — cos/3)°/"' 

X P„_i/2(cosha<)Q„_i/2(coshQ;>). (25) 
Now, using cosn(/3 — /3') = cos n/3 cos n/3' + sinn/3sinn/3', we obtain 



^ 6Mssmh _ cos^' ^ eoPn_i/2(cosha<)Q„_i/2(cosha>) 



27r 

n=0 

/■^'^ sin;3cosn;3d/3 ^, Z"^'' sin/3sinn^d^ 

cos rtp / , + smnp 



(cosh a — cos /3)^/^ Jq (cosh a — cos /3)^/^ 



(26) 



The first integral vanishes identically (integrand is odd under ji — > — /3), whereas the second one 
gives (Ref. [23], p. 961): 

^ 2- sin Sinn/3 ^ 2n cosn/3 ^ _ Jl^gi (27) 

(cosh ao — cos /3)^/^ 3 Jo (cosh ao — cos /S)-^/^ 3sinhao "^^/^ 

Now, once we are dealing with a toroidal shell, so that ol — a — a< = q;> = ao, we finally 
obtain: 



4 V26Ms sinh ao , flM/2V^ • o 

•Pjr = (cosh QfQ — cos p) ' > ne„sm np 

37r ^-^ 



X -Pn-i/2(coshQ;o)(3„_i/2(coshao)Q;^_i/2(coshao). (28) 
Now, substitution of (1^51) in gives: 

^ 4^/2^oA^i&^ , 3 

-^mag = ^ Smh ao 2^ "-enf»i-l/2(yn-l/2 

n=0 

1 /■^'^ sin/3sinn^d/3 
^"-1/2 7^ (coshao-cos/3)5/2' ^^^^ 

which after integration by parts immediactly yields Eq. ([5]), as desired. 



15 



C.2 Evaluation of Eqs. ([9]) and (fTOb 

For the onion-type patterns, we have: 

cosh a cos /? — 1 



M ■n = Ms 



mr{ip), 



and 



cosh a — cos /3 
/ cosh a — cos /3 \ dm^{ip) 



\ b sinh a 

Now, the substitution of ([301)- ([311) in ©, we obtain: 



dip 



(30) 



(31) 



$o = ^(cosha--cos/30^/^Ef:(-l^^^" w! , U ^^-i 



47r2 



fc=0 n=0 



T(n + A:- 



/2 



x<^ sinhaoQ„_i/2 / d/3 



27r 



J (cosh ao cos (3—1) cos n{/3 — /?') 
(coshao — cos/3)^/2 



mr{ip) COS fc(iy9 — ip')dip 



Qn-i/2icosha) /-s- cosn(/3-/3') f^^ dm^{^) 
-da 



sinh a 



(cosh a — cos /?)3/2 



■ cos k{Lp — ip') dip \ , (32) 



with (coshao) = P/f, and similarly for Q. Using the usual identity cos7i(/3 — /?') = 
cosn/3cos7z/3' + sinn/3sinn/3', and integration by parts, we get: 



2V2bA'L 



= J^"^"'' (cosh a' - cos ,9')^^^^ y^(-l)''e„£fc cos n/3 



^,r(n-fc + i) 



37r2 



r(n + fc + i)^"-i/2 



fc=0 n=0 



(33) 



where 



cosh ao 



sinhao "^"-1/2) sinh 



ao 



Jki'p') ~ / mr(iy9) COS fc((y9— iy9')(i(^ 



3 /-"^ Q^_i/2(cosha)Qi_-^,2(cosha) 



sinhf 



da, 



cosk(p—ip )dp. 



dip 



After a long and tedious algebraic manipulation, substitution of expression ([33]) in Eq.© yields 
the desired Eq. ([9]), from what Eq. (fTO|) is readily obtained by setting q — 1. 
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